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Abstract 

We analyze the heavy quark bound state spectrum using an order-dependent 
conformal mapping to re-sum the perturbative expansion for current corre- 
lators. The procedure consists of two main steps. Firstly, the Borel plane 
structure of the truncated perturbative expression is modified to ensure con- 
sistency with the operator product expansion. This is analogous to a resum- 
mation of infrared renormalon chains. Secondly, this perturbative expansion 
is conformally mapped to a new series with improved convergence properties. 
This approach may be shown to induce power corrections consistent with ex- 
isting condensates, and the resulting expansion may be ordered in powers of 
an infrared-analytic effective coupling. The technique is then applied to cc 
and bb sum rules without any explicit introduction of vacuum condensate pa- 
rameters. Ground state masses for the vector, axial-vector and A' channels 
are well reproduced, while results for the scalar-pseudoscalar mass splitting 
are less impressive. 
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I. INTRODUCTION 



The Wilson operator product expansion, within the framework of QCD sum rules [1]-|3| 
(see also || for a review), constitutes the standard continuum approach for the system- 
atic treatment of nonperturbative corrections to perturbation theory in QCD, and study of 
the bound state resonances. Taking asymptotic freedom and perturbative QCD as a start- 
ing point, this formalism extends the validity of the operator product expansion (OPE) to 
the scale at which resonances appear in the spectrum as a manifestation of confinement. 
This generalised operator product expansion for the correlators of various electromagnetic 
currents includes, in addition to the perturbative contribution, non-perturbative power cor- 
rections associated with vacuum condensates of the fields. 

For a given current j r , we have the correlator^], 

T^...IT r (Q 2 ) = i I d 4 xe^ (0\T{jI..(x)jI.(0)}\0) , (1) 

where T^... is a tensor which depends on the Lorentz structure of the particular current 
j r . The validity of the OPE at the resonance scale requires the assumption that one may 
separate the short distance (perturbative) and long distance (condensate) contributions in 
a consistent manner. The expansion then has the form, 

i J d^T(jl...(x)f u ...(0)) = Cl ert (Q*)l + £ ^P(?2„, (2) 



n>2 



where Cp ert (Q 2 ) is the perturbative contribution, and Oi n are 2n-dimensional gauge invari- 
ant operators with C^y the corresponding perturbatively calculable Wilson coefficients. 

While the perturbative contribution is calculable in a consistent manner in the short 
distance regime, provided one imposes an infrared (IR) cutoff, the fact that the split is 
not necessarily unique or unambiguous is signalled by the perturbative large order diver- 
gences associated with infrared renormalons [p|-[T0[|. Indeed, it is well known that the Borel 



resummation ambiguity of the perturbative series for correlators associated with inclusive 



processes has precisely the form, if not the magnitude ||11|| , of the power corrections de- 
termined by the non-perturbative condensate contribution to the OPE. The perturbative 
ambiguity may then be cancelled by a similar ambiguity in the condensates, rendering the 
combination well-defined. Consequently, estimating the numerical value of condensates re- 
quires detailed knowledge of the perturbative ambiguity itself, hence the recent interest in 
renormalon asymptotics (see e.g. |fL2"l - |24}l and references therein). 

Our aim in this paper is to consider the calculation of the mass of heavy quark bound 
states within the sum rules formalism via an alternative approach which is to try to resum 
the asymptotic perturbative series into a convergent form. There are two clear problems 
that arise in this endeavour which need to be overcome: (1) The Borel non-summability of 
the perturbative series, due to the presence of IR renormalons, implies that such a resum- 
mation technique must be powerful enough to induce power corrections to counteract the 



1 We define Q 2 = —q 2 , with q the current momentum. 
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fundamental ambiguity in the asymptotic perturbative series; (2) The Borel resummation 
ambiguity, while having the correct functional form to hint at the structure of the nonper- 
turbative corrections required, cannot itself provide the normalisation of these terms from 
within perturbation theory, as the condensates have a non-perturbative origin associated 
with confinement to which perturbation theory itself is ostensibly blind. 

Nevertheless, with regard to point (1), order-dependent mapping techniques, based on 
the resummation of a perturbative series in the coupling g, G(g), via a conformal mapping 



to a new series in a, G(g(a)), p5| , p6| are powerful enough to resum Borel non-summable 
series into convergent form. An example is the quantum mechanical double-well potential, 
and we illustrate the structure of the induced corrections in a 0q model in Appendix A to 
indicate their connection with field theoretic power corrections. 

The nonperturbative input stressed in point (2) can, however, only arise in such an ap- 
proach in a summation to all orders, and thus cannot obviously be obtained with knowledge 
of only a finite truncation of the perturbative series, which is the practical reality in field 
theoretic situations. Furthermore, rigorous proofs of convergence require detailed knowl- 
edge of the complex- analytic structure of G(g), which is again generally lacking in 4D field 
theories. Thus, in this case one needs extra information to try and "optimise convergence" 
of the conformally mapped series. This is possible since in Eq. @ there is a free positive 
parameter C (a is usually fixed via constraints on the mapping) which may be fixed order 
by order to improve convergence. Generically, a necessary condition is that C must scale 
with the order of the expansion in a specific manner. In QCD, with our knowledge limited 
to only the lowest order terms, a natural way to achieve this is to use some infrared data 
related to confinement, as this non-perturbative input allows us to estimate the appropriate 
scaling of C for the resummed series. Note, however, that our implicit parametrization of 
the nonperturbative corrections which may be relevant is directly related to the structure 
of the conformal mapping, and not just the constant C itself. 

An order- dependent mapping approach of the kind discussed above has been developed 
by Solovtsov p7| , p8[] (see also p9|-p4| for applications) which we briefly review in Section 2, 
indicating how, at any finite order, the conformally mapped series may be re-ordered into a 
form structurally equivalent to perturbation theory, but with an infrared analytic effective 
coupling which alters the IR behaviour of correlation functions. 

For the application of this technique in situations where the numerical value of condensate 
related power corrections is significant, and only one- and two-loop perturbative coefficients 
are known, as in the sum-rules case, we also require a modification of the perturbative result 
to ensure the correct all-order resummation ambiguity, i.e. the correct position of the first 
IR renormalon pole. This then ensures consistency with the OPE prior to resummation 
in the sense that the perturbative ambiguity may be consistently cancelled by ambiguities 
in the OPE condensates. Or in the present context, that any induced power corrections 
compensating the perturbative ambiguity can be associated with OPE condensates. An 
important point is that this modification is only necessary to ensure the correct momentum 
dependence of the ambiguity, and thus a full resummation is not required. The normalisation 
is provided by the nonperturbative parameter C, and not the renormalon residue. 



a 



(3) 
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We illustrate this idea in the massless case in Section 2, which was also discussed in 
34], while in Section 3 we apply the technique in the determination of the masses of the 
cc and bb bound state families. In this case we also require a resummation of threshold 
Coulomb singularities and the particular resummation procedure adopted here is discussed 
in Appendix B. We then show that one may closely approximate the experimental estimates 
in all channels other than the scalar and pseudoscalar. This suggests additional sources of 
power-behaved contributions in those channels to which this technique is insensitive. In 
Section 4, we discuss our results and consider the possible interpretations. 



II. AN ORDER-DEPENDENT MAPPING IN QCD 

In this section we briefly review how an order-dependent mapping of Green functions may 
be obtained via direct manipulation of the functional integral in QCD using the approach 
of Solovtsov I^^SJ (see also P§|-|3"4j]), and then proceed to discuss the developments which 
facilitate later application to heavy quark sum rules. 

At a formal level, the mapping may be implemented along the ideas of Seznec and 



Zinn- Justin |25[ via a direct resummation of a perturbative series G(....) = I]„c n A n 
expansion in a new parameter a, G = G(X(a)), related to A = a s /(47r), by the conformal 
mapping, 

A = - - a " (4) 

(4vr) 2 C(l-a) 3 ' [) 

where C is a positive constant. Note that this is particular case of the conformal mapping 
(H) with a = 3/2, and one observes that < a < 1 for all values of the gauge coupling. 

However, in the case of QCD, to ensure gauge invariance, and to allow implementation 
of the renormalization group, it is convenient to have an explicit realisation of this mapping. 
Such a realisation has been developed by Solovtsov |27|j28| , and becomes possible with the 
introduction of an auxiliary field used to split up the quartic gauge field interaction 
term, which we denote ^(A)!, 



exp(tg 2 S 4 (A)) = J [d X ] exp 1 J dx X % [x^ a + t^=f abc A^ b A 



(5) 



Use of this auxiliary field reduces all the interaction terms in the QCD action to Yukawa 
form, and enables a new split between free and interaction parts, S x = S + Sj, parametrised 
by the variables £ and £, where 

S = C'ISiA, X ) + Siffl + S 2 (c)} + r 1 S 2 (x), (6) 

and 



2 We note that the theory is still consistent off-shell. One may readily verify that if x has a gauge 
transformation consistent with its on-shell constraint, x ~ 9fAA, then the full effective action still 
satisfies the functional Slavnov- Taylor identities. 
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Sj = gS 3 (A, ^ c) - (C 1 - 1)[S(A, X ) + + S,(c)] - (r 1 - l)^*)- (7) 

We use a compact notation for the standard kinetic and Yukawa interaction terms of the 
gauge (A), quark (ip), and ghost (c) fields, and there are also terms fixing the covariant «g 
gauge. if>2(x) is a mass term for the x ne ld introduced in (|), while S(A,x) is the gauge 



propagator in the x background (see |28|). A Green function for an even number of fields 
then takes the form, 



<?(....) = ^ fa^^cjlElG-^W^, (8) 

Z J n=0 n[ 



where [d(A, if), c)] denotes the conventional functional integral over gauge, quark, and ghost 
fields, and Z is the partition function. 

As discussed in [pT[p8"| , further expansion of Sj, and a rescaling of the fields, allows x 



to integrated out, restoring the standard free action of QCD in the exponential and the 
correct relationship between the 3-point and 4-point gauge couplings, provided that £ = ( 3 . 
However, on inspection of the subsequent series, one observes that it may be ordered in 
terms of a new parameter, a = 1 — (, provided one identifies a with the gauge coupling 
precisely via the conformal mapping @. 

Expansion of G as a power series in a, up to 0(a 5 ), then results in the expression 

G (5) (-) = <*>(..) + ^ 2 (..) + 3^ 2 (..) + ^[QCg 2 (..) + g,{..)\ 

+ ^[10^ 2 (..) + 6^ 4 (..)] + O(a 6 ), (9) 
where g2n{--) are the corresponding terms in a standard perturbative series with coefficient 



A n [pq| . Re-ordering this series in terms of the perturbative factors g2n(--) we find 



G( 5 )(..) = go (..) + £(1 + 3a + 6a 2 + 10a 3 



■■■)92{: 



a 4 



+— (l + 6a + ---)^(..)] + --- (10) 
= 0o(-O + A e/ /i 5 ^ 2 (-) + (\ ef f) 2 \M..) + ■■■, (11) 

where we have introduced the effective coupling, A e //| n , which is the expansion of (|) to 
0(a n ). Therefore, at a given finite order, the conformal mapping may be realised by a direct 
replacement of the coupling A in the perturbative series with A e // defined to the appropriate 
order in a. We shall henceforth work at 0(a 3 ), and thus 

A e// = £(l + 3a). (12) 

While the equation ([|) implies a formally nonlinear relationship between the power series 
in A and a, for the truncated series truly non-perturbative input arises via the ability to fix 
the variational parameter C at each order. One thus obtains a sequence of approximants 
for the quantity G under study {Gi(0(a), C{), G 2 (0(a 2 ), C 2 ), . . . , Gn(0(ci n ), Cn)}, and the 
choice of C is dictated by the aim of optimising the convergence of this sequence. 
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In the present case the normalisation of the result depends on the value of condensates 
related to confinement effects which are invisible at any finite order of perturbation theory. 
The structure of these corrections is implied by the large order ambiguity, but this does not 
assist in providing the correct normalisation. Therefore we fix C at each order by fitting a 
particular low-energy quantity, which ensures the correct normalisation at that scale. The 
universality of this choice can then by checked by applying the same method to different 
processes, and the relative magnitude of successive terms in the series can still be monitored, 
even though a formal proof of convergence is lacking^. 

Given this philosophy, the approach which has been used (with previous success PP|-P^[) 
is to ensure that the running coupling has a singular infrared asymptotic behaviour consis- 



tent with the linear part of the static quark potential, Vu n (r) = ar. As discussed in p8 
this corresponds to the requirement that the /3-function obey /3(A) — > —A for large A. 
Enforcement of this constraint requires us to consider the renormalisation group evolution 
of the parameters. A convenient check on the gauge invariance of the expansion is provided 
by calculating the charge renormalisation constant with an arbitrary covariant gauge and 
ensuring that dependence on the gauge parameter drops out at each order in a. This is 
indeed the case p8fl , and to 0(a 3 ) the /3-function is given by 



gM -^ (2 + a Hl-a)3 < 2 + 9a >- < 13 > 

where b is the first coefficient of the perturbative /3-function. The resulting RG equation 
has the form 



2& 0l Q 2 



/(o) = /(ao)+^ln^, (14) 



where at this order 

2 6 , 18 1 624 , , . 5184 , , 9 

/ (a ) = 48 In a - - + — In 1 - a + — - In 1 + - 

a z a 11 1 — a 121 121 2 



The parameter a and the momentum Q in Eq. (14) are defined at a normalisation point for 
the effective coupling constant (|i~2"l). Explicitly, using phenomenological data for the string 
tension a, the above analysis of the /3-function has been shown to lead to C = 4.1 [[28[] at 
this order. The running coupling is then given by (|T|), with the running parameter a(Q) 
determined implicitly by ([14]). As noted earlier, a(Q 2 ) is finite at all scales and thus A e // 



does not exhibit a Landau pole [28 



Having fixed parameters in the infrared, the first consistency check is that we recover 
the perturbative expression for the running coupling at large energy scales. Indeed this is 
the case. If, using ([P2|), ([14]), and (|i~5|) , we consider the regime Q 2 ^> Ql = A.q CD where 
A ~ a 2 /C 1, we recover the one-loop perturbative result, \(Q 2 ) — * 1/ {bo\n.{Q 2 / K 2 q CD )). 



3 Note that such series have been proven to converge in simpler models p5|-|37[] for which the 
corresponding perturbative series is at best asymptotic, and also as discussed in Appendix A in 
some situations where the perturbative series is non-Borel-summable [26]. 
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The IR running of the coupling is, however, significantly modified [ 28| . The running quark 
mass [27], given by 



m{Q2) - ™° (18) 10 • (16) 

where mo = m(Ql) and 70 = 4 is the first coefficient of the anomalous dimension, similarly 
coincides with perturbation theory at large scales but has a modified IR behaviour. 

We now consider the application of this approach to the current correlators introduced 



in section 1. For illustration it is convenient to consider first massless correlators [51], and 
in particular the Adler D-function, D(Q 2 ) = — Q 2 dU(Q 2 )/dQ 2 . These arguments will then 
be extended to the massive case in Section 3. 

We begin with the standard spectral representation, 

poo 1 

D(Q 2 ,X) = Q 2 j o ^s^^R^X), (17) 

where R{s) = ImU(s + ie)/ir is given in a convenient normalisation at 0(A) by 

R{s) = 1 + 4A. (18) 

Conventional RG improvement of the D-function would break its analytic properties by in- 
troducing the perturbative Landau pole at Q 2 = Aq CD . However, this may be circumvented 
within the spectral representation using RG improvement of the integrand ||38|| , with the 
knowledge that R{s) obeys the same homogeneous RG equation as D. This corresponds to 
choosing a more general solution of the RG equation for D with the same UV asymptotics. 
However, the integral in ([17]) is then undefined, as it now runs over the Landau pole. The 
structure of the virtuality distribution in flTTD implies that there is an infrared renormalon 
ambiguity of 0(Aq CD /Q 2 ). However, there is no corresponding condensate operator in the 
OPE, and thus this first IR pole is expected to be absent, at least within the current per- 
turbative analysis. An all-order resummation in the large &o limit would indeed recover the 
expected branch structure, removing the first IR renormalon pole. 

We shall make the assumption that the first IR pole should be absent from the point 
of view of perturbative asymptotics, to ensure consistency with the OPEQ However, as 
discussed in Section 1, a full resummation is unnecessary in the present context, as the 
normalisation is fixed elsewhere. Thus we can use a simple trick to remove the first IR pole 
and ensure that the large order asymptotics are consistent with the OPE. A convenient way 
to achieve this is to perform an integration by parts in fli~7|) [Q, obtaining 



4 Note that this excludes the possibility that such a correction could arise as an exponentially 
suppressed correction to the coefficient of the identity operator in the OPE. 
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FIG. 1. The virtuality distribution functions tw(t) taken from Ref. [^3| (solid line) and the 



function (21) multiplied by a factor of r (dashed line) versus lnr. 



At the order to which we are working the second term vanishes, and we find 

D(t,X) = 1 + 4 / d™(r)A(/z 2 ), 
Jo 

where t = Q 2 / n 2 . The weight function, given by 

W (r) 



2r 



1 + r 



(20) 



(21) 



describes the distribution of virtuality which we shall comment on shortly. Performing RG 
improvement on the integrand, using for consistency the one-loop /3-function f3(X) = —boX 2 , 
we obtain 



/•CO 

D(t,X) = 1 + 4 / cWr)- 
Jo ' 1 



A 



Xb n In r 



(22) 



where A = A(t).n The modified structure of the virtuality distribution, which at low scales 
now has a linear dependence on r, u(t) ~ 2r + 0(r 2 ), is now consistent with an ambiguity 
associated with the first IR renormalon pole of 0(Aq CD /Q a ) which can be consistently 
cancelled by an associated ambiguity in the gluon condensate in the OPE. 



We observe that the virtuality distribution function fl2lD coincides with the function used 
in and remarkably is numerically very close to that obtained in P3| (see Figure 1) which, 
in contrast to the present naive modification, corresponds to an all orders resummation of 
renormalon contributions in the large 6q limit. 



5 We denote perturbative running parameters in the MS scheme with a bar, as opposed to the 
unbarred running parameters in the variational series. 



S 



This connection with renormalons is illustrated more clearly by performing a formal 
Borel transformation on (0), whereby we obtain 



(23) 



with 

B(6) = r(l+66 )r(2-66 ). (24) 

This Borel function exhibits the correct infrared and ultraviolet renormalon poles for the D- 
function, although not the full branch structure flO ]. Nonetheless, as mentioned above, the 
crucial point here in using the partial integration in (|2T|) is to obtain the correct positioning 
of the poles and thus an ambiguity consistent with the OPE. 

Having massaged the perturbative series into a form whose large order divergence may 
be consistently cancelled by ambiguities in the OPE condensates we may now perform the 
conformal mapping, as discussed earlier, by replacing the perturbative running coupling with 
X e ff (|l|). In the modified spectral representation, this leads to the replacement of (|22|) by 

D(t, A) = 1 + 4 / drw(r)Ae//(tr), (25) 

J 

where, to 0(a 3 ), A e // is given by (|l|), with the running expansion parameter a = a(Q 2 ) 
determined via fll4"|). 

Within a field theory such as QCD it is not possible to directly assess the convergence 
of this series, and therefore to determine whether indeed OPE type power corrections of 
the appropriate magnitude are induced when compared to the initial perturbative series^. 
Although analysis of a 0q model |2(J (see Appendix A) indicates that the mapping is able 
to resum the appropriate corrections, within QCD it is not always possible to assess even 
the magnitude of higher order terms. Therefore we propose instead to test the technique by 
analysing the massive case, and comparing the results directly with experiment. 



III. MOMENTS FOR HEAVY QQ BOUND STATES 

We shall make use of the standard sum rules approach for the determination of the lowest 
mass resonances |lpf| (see also [§]] for a review), wherein one assumes the validity of the 
narrow resonance approximation for the lowest mass contribution to the imaginary part of 
the two-point current correlator. The other side of the sum rule is conventionally deter- 
mined directly from QCD up to parameters associated with vacuum condensate operators. 
However, in the present context, as elaborated earlier, there is no distinct split between 
perturbative and nonperturbative contributions. The non-perturbative parameter C does 
not have any direct connection with vacuum condensates. 



6 Note that a 1/Q 2 correction may also be induced in the running coupling by removal of the 
Landau pole |l8|]. However, this has a purely perturbative short distance origin. For recent work 
see Hill. 
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We now repeat the analysis of the preceding section in the case of massive correlators. 
To 0(A) in perturbation theory we have 



Irnlfh) = — 
Air 



uP(s)+4XU^\s) , (26) 



where T denotes the current in question, and s = q 2 . The one- (IT )) and two-loop (IlW) 
components have been obtained previously and are enumerated in Appendix B to fix our 
notation. The first convergent moment is defined by j§] 



du r (Q 2 ) 

' J - V) ^ ' dQ 2 



ds 



TT J im 2 ( S + Q 2 ) 2 



/mli r (s), (27) 



which is also the Adler D-function up to a factor of Q 2 . Introducing o = s — Am 2 , which we 
associate with the virtuality, and the quantity u 2 = o jio + 4m 2 ), we have 

An additional problem encountered in studying moments for heavy qq systems is the 
dominant Coulomb interaction, which leads to the perturbative series having an effective 
expansion parameter given by X/u. In the conventional approach, higher order moments 
are dominated by the nonrelativistic region of small u and thus the perturbative expansion 
breaks down due to the size of the expansion parameter ^T|. Since our approach still makes 
use of these perturbative coefficients a resummation is required in order to consider the 
moments for large n. 

When, in analogy with the massless case, we perform an integration by parts to remove 
the first IR renormalon pole, a particular resummation which places higher order corrections 
inside the virtuality distribution becomes very natural. This resummation may be performed 
exactly, and in Appendix C we describe its implementation in a simplified model illustrating 
its effectiveness in resumming Coulomb singularities, and its connection with the familiar 
Sommerfeld-Sakharov factor p2 |. 

Implementation of the resummation in the present case results in the following expression: 

m wq) = ^/ d ° {Q2+W{a)r+Nr > (29) 



where to 0(A) 



W (a) = (a + Am 2 ) ( 1 - 4A^#- ) , (30) 




and 



Mu) = (l-u 2 r»r£ du ^ ' U?(u). (31) 
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Carrying out renormalisation group improvement of the integrandQ, the solution of the 
RG equation corresponds to ( p9|) with A replaced by A(cr) and m replaced by m(a), where in 



the MS scheme the arguments of these running parameters are scaled by k\ = exp(— 5/3) 

The conformal mapping is then realised via the replacement of perturbative running 
parameters with the non-perturbative running parameters in (|l^) and (|16D , and in place of 
(P0|) we have 



W(a) = (a + Am\ka)) ( 1 - 4X eff (ka)^0- ) , (32) 

ll r [u) , 



where now u 2 = cr/(cr + Am 2 {ka)). 

Finally, we may obtain the power moments for arbitrary n, which are given by 



°°, (a + 4m 2 {ka)) N ^IlP' 



Ait 2 Jo (Q 2 + W(a)) 1+n 

As in the standard sum rules approach, the mass of the first resonance is then obtained by 
considering ratios of the moments, = M^i/M^, for large n. 

Using the perturbative formulae for the imaginary parts of the two-point correlators 
listed in Appendix B the integral for the moments can be evaluated numerically for each 
current. There are only three parameters in this calculation: a reference value for the QCD 
coupling ao, taken at the r mass scale; the quark mass m c (or m^); and the variational 
parameter C, constrained to be near 4.1 via data from the quark- ant iquark potential. We 
have previously presented an analysis of r-decay using this technique |54[] which resulted in 



a particular normalisation of the coupling at the r scale, ao(M r ). For consistency we use 
the value extracted at 0(a 3 ), «o = 0.379, which corresponds to the order of the perturbative 
expressions with which we are working. 

An important aspect of the modification of the virtuality distribution (Fig. 1) invoked to 
ensure consistency with the OPE, is that the peak is shifted to higher scales. Consequently, 
implementation of the Coulomb resummation produces a saddle point for the moments which 
dominates the expression for large n (see Appendix C). The existence of the saddle point is a 
result of the function W (a) attaining a minimum, at say o = b. This scale a then dominates 
the expression for the moment as n becomes large. Analysis of the various channels indicates 
that this minimum occurs in the range a ~ 4 — 6 GeV 2 for cc, and a ~ 20 — 30 GeV 2 for bb. 
The fact that the saddle point occurs well above Aqcd provides additional justification for 
the validity of the perturbative expressions with which we are working. It is then clear that 
for large n, the moment ratio tends asymptotically to its saddle point approximation, i.e. 

R n (Q 2 ) n ^Q 2 + W(a). (34) 



7 Note that for currents with a non-zero anomalous dimension the modification takes the form of 
an overall factor in the integral which will not contribute to the moment ratios for large n, and 
thus not to the asymptotic estimates for the bound state masses. For this reason we may ignore 
any anomalous dimensions. 
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The corresponding ratio, R^iQ 2 ), which arises from assuming the validity of the narrow 
resonance hypothesis, has the form Q 2 + M R for large n, where Mr is the mass of the first 
resonance in the relevant channel. Thus we obtain the prediction 



M R = s/W(a), (35) 

and as a consequence there is no requirement to fix Q 2 explicitly, for example at Q 2 = 0, in 
order to determine the mass of the resonance. 

Results of the numerical calculations for the full moment ratios, and the asymptotic 



estimates yW(a), for all currents under consideration and for both cc and bb bound states 
are shown in Figures The quark mass parameters giving the optimal fit for all channels 
in each family are m c = 1.51 GeV, and m& = 4.72 GeV, which agree well with experimental 
constraints ||44|| . As a consequence of the resummation, we observe that in all cases the 



moment ratios are stable as n becomes large, in contrast to the conventional results obtained 
by explicitly adding power corrections to a truncated perturbative series for condensates up 
to dimension six or eight. For clarity a summary of the results presented in the figures is given 
in Tables 1 and 2, where the asymptotic estimates for the moment ratios are compared with 
the experimental results [M\ where available. For the vector, axial-vector, and A' channels 



the magnitudes of the bound state masses, and also the inter-channel splitting, are very 
well reproduced. For the scalar and pseudo-scalar channels the magnitudes are reasonably 
well approximated but the splitting between these states is not well described. Possible 
interpretations for the results obtained will be presented in the next section. 
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n <j(GeV=) 

FIG. 2. On the left we plot ratios of the vector (solid curve), axial vector (dashed curve) and 
A' (dot-dashed curve) moments for cc bound states, and on the right we plot y/W(a) versus a for 
the same currents, the minimum being the asymptotic limit of the moment ratios for large n. For 
comparison, the straight lines are, in order of decreasing mass, the corresponding experimental cc 
A', axial- vector, and vector bound state masses p4| . 
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a = 0.379 

















Parameters: C = 4.1 

m =1.51 
a„=0.379 



5 

u(GeV=) 



FIG. 3. On the left we plot ratios of the pseudo-scalar (solid curve), and scalar (dashed curve) 
moments for cc bound states, and on the right we plot \/W(a) versus a for the same currents, the 
minimum being the asymptotic limit of the moment ratios for large n. For comparison, the straight 
lines are, in order of decreasing mass, the corresponding experimental cc scalar and pseudo-scalar 
bound state masses 44]. 
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FIG. 4. On the left we plot ratios of the vector (solid curve), axial vector (dashed curve) and 
A' (dot-dashed curve) moments for bb bound states, and on the right we plot y/W(a) versus a for 
the same currents, the minimum being the asymptotic limit of the moment ratios for large n. For 
comparison, the straight lines are, in order of decreasing mass, the corresponding experimental bb 
bound state masses for the axial- vector and vector channels ^4j . 




FIG. 5. On the left we plot ratios of the pseudo-scalar (solid curve), and scalar (dashed curve) 
moments for bb bound states, and on the right we plot y/W(a) versus a for the same currents, 
the minimum being the asymptotic limit of the moment ratios for large n. For comparison, the 
straight line is the corresponding experimental scalar bb bound state mass [44]. 
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TABLES 



current 


M expt (GeV) 


y/W(a) (GeV) 


a (GeV 2 ) 


3V 


3.10 


3.06 


4.6 


3A 


3.51 


3.52 


4.8 


3A> 


3.53 


3.56 


4.8 


3S 


3.42 


3.27 


6.0 


3P 


2.98 


3.19 


3.9 



TABLE I. Summary of the estimates obtained for the cc bound state masses, compared to the 



experimental values [44]. Note that a denotes the virtuality scale at which the asymptotic result 
was obtained. 



current 


M expt (GeV) 


y/W{a) (GeV) 


a (GeV 2 ) 


3V 


9.46 


9.43 


27 


3A 


9.89 


10.02 


26 


3A> 




10.07 


26 


is 


9.86 


9.80 


30 


jp 




9.56 


24 



TABLE II. Summary of the estimates obtained for the bb bound state masses, compared to the 



experimental values [44], with a again denoting the virtuality scale at which the asymptotic result 
was obtained. 
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IV. DISCUSSION 



In this paper we have presented an explicit application of the order- dependent map- 
ping approach to the study of the heavy quark QCD resonance spectrum. While it proved 
convenient to study and modify the perturbative asymptotics to ensure consistency with am- 
biguities in the OPE power corrections, it is clear that subsequent to the conformal mapping 
any obvious split between perturbative (short distance) and nonperturbative (long-distance) 
contributions is lost. Thus, other than by fitting data, it does not appear possible to extract 
values for the condensates themselves or even to determine the dimension of the dominant 
contribution, in contrast to related approaches such as that of Dokshitzer, Marchesini, and 
Webber ||45|| . This requires a clear and well-defined split between the perturbative and 
nonperturbative contributions, a distinction which does not emerge naturally within the 
framework of order-dependent mappings. 

Nevertheless, at a purely calculational level, applying this formalism to the bound state 
spectrum, without the need for explicit introduction of condensate parameters, we found 
quite remarkable agreement with experimental results for the vector, axial-vector, and A' 
channels in both cc and bb families. Indeed if we use a more accurate extraction of the 
coupling at the r scale at 0(a 5 ) |34| the results are even more impressive for these channels, 
although for consistency one should then include the next order perturbative coefficients. 

However, the results for the scalar and pseudo-scalar channels, while again having ap- 
proximately correct magnitudes, were unsatisfactory in terms of the relative splitting be- 
tween the states. This was not improved by a more accurate extraction of the coupling. It 
may be that higher order corrections within this approach could resolve this discrepancy. 
However, considering the success in the case of the vectorial channels, it appears likely 
that there are contributions to the scalar and pseudo-scalar channels which are missed by 
the variational approach, e.g. matrix elements not associated with vacuum condensates. 
Whether this is due to the particular averaging over non-perturbative contributions implicit 
here, or more closely related to use of the static quark potential model for normalisation, is 
currently under investigation. 
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APPENDIX A: 



Induced Power Corrections in a (p% =0 Model 

In this appendix we use a simple zero dimensional model discussed by Guida, Konishi, 
and Suzuki |2B| to illustrate concretely the power of order-dependent mappings in gener- 
ating power corrections to a perturbative series, and consequently producing a convergent 
expansion in a situation where the corresponding perturbative series is non-Borel summable. 

Consider the zero dimensional 4>% =0 model with partition function 

/OO OA 
d^e-f 2 -^, (Al) 

and analytically continue the result for g gC. With the change of variables z = y/2g<p 2 , this 
is represented in terms of a parabolic cylinder function D v (a) |f46f , 



An N th order approximant to this partition function in the order-dependent mapping 
may be written in the form 

Z N (a) = {l-a) l l 2 ip N {a,C N ), (A3) 

where ipN is the A </l -order truncation of the Taylor series for a function ip. The order- 
dependent mapping between the coupling g and the parameter a is given by 

9 = cfrbr (A4) 

where this relation holds for complex g. Restricted to the positive real axis this N th order 
approximant may be written explicitly as 

Z N (a) =ff^(l- af' 2 a k , F ( n + k + 1 / 2 ) (A5 ) 

where u n are the perturbative expansion coefficients 
Convergence of the sequence of approximants, 

{Z N (a,C N )} N -^?Z(g) (A6) 

has been proven in [^6| for the analytically continued model, and the details will not be 
recalled here. The importance of this result in the present context for investigating induced 
power corrections follows from the structure of the asymptotic series for Z(g) for \g\ 1. 
This expansion has the form 

^ £ r(2n + 1/2) ((-!)" + rmV^V, (AT) 

n=0 n - 
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with m = ±1 for =p57r/2 < Arg(g) < =F7r/2 and m = otherwise. If we consider the 
structure of this expansion for g = e l7T g R with g R real and positive, then the first term is the 
standard perturbative series about g R = 0, which at large orders has coefficients exhibiting 
the factorial growth ( J e / '2tcT \n)2 2n ) associated with a Borel non-summable series. The 



second term can be interpreted as a complex nonperturbative "power correction" accounting 



for the imaginary part introduced in resolving the perturbative ambiguity (c.f. [15]). 

This may be seen more explicitly by performing a Borel transform on the perturbative 
series for Z{e m g R ). From 

oo n 

Z pert (e™g R ) = ]T ^fr(2n + 1/2), (A8) 

representing the T-function as an integral and transforming variables, we have formally 

IT r A 
9r Jo 



Z per \e i7r g R ) = J— f dbe- b/9R B(b), (A9) 
V Or Jo 

where 



B{b) = (A10) 



exhibiting branch points at b = and 6 = 1/4. Strictly, the transformation of variables we 
have performed is only valid for b < 1/4 and some regularization is required to integrate 
beyond A = 1/4. Nonetheless, the position of the branch point at 1/4 can also be inferred 
directly from the structure of the large order coefficients Je/2TiT(n)2 2n . 



The branch point at b = may be removed by renormalisation, while the ambiguity 
associated with regularising the integral at the branch point b = 1/4 is given by the residue 
as (c.f. the second term in (]A7])) 



Res{Zr:t /A ) = J-e-- R . (All) 
V 9r 

For illustration, we may symbolically represent this ambiguity as it would appear in a field 
theory with 

,J « m = wgyvy (A12) 

where /3(g R ) = —bog R + ■ • -. The ambiguity then takes the form 

exp (-iH!f 4 - <Ai3) 

which we may loosely interpret as an OPE type "power correction" . 

Thus we can conclude that the weak coupling expansion is only consistent when one 
includes a power correction of this form to the perturbative series. More generally the 
regularisation of the Borel singularity will also result in an imaginary contribution (e.g. by 
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passing the contour above or below the branch point) counteracted by the imaginary part 
in (|A7|). 

Therefore at least in this rather formal example, the proven convergence of the order 
dependent mapping ensures that it can generate a power correction of this form. Note 
that it is likely that one can perturb this result away from the branch cut by considering 
Z^e^^'gji). The imaginary part of the weak coupling expansion still exists for any e < 7r/2, 
so one might expect that a power correction is still required in this case. 
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APPENDIX B: 



Perturbative Wilson Coefficients 



For completeness, and also to fix our notation, in this appendix we list the perturbative 
Wilson coefficients for the relevant correlators. The currents considered are listed below in 
the format jr = ...(J PC ): 

• scalar : j s = ipi4>j (0 ++ ) 

• pseudo— scalar : jp = iipfjgipj (0 h ) 

• vector : j v = ipa^j 

• axial vector : j A = {q^/q 2 - g^)4>flul^j (1 ++ ) 

• A' : j A , = ^<9 M 7 5 V>j (l +_ ) 



Using the notation of we introduce the following generic components for the correla- 



tors, 



A(u) = (1 + u 2 ) 



7T 



h In 

6 L — U 



l+u, l+u 



In 



+ 21 



A'(u) 



[l + u 1 



-21 



21 



u 



Al(u) + l{u 2 ) 



1 - u 



u 



+ 3u In ■ 



+ 21 



u 



hi 



l + u 



— uluu 



u 



-21 



u - 1 
u + 1 



1 — u l + u 1 — u 
—3 In In h 2 In In u 



l + u 



l + u 



where 



Ux) = - f X dt-ln(l -t), 
Jo t 



(Bl) 



(B2) 



(B3) 



is the Spence function. Then following the notation of Section 3, the relevant formulae for 
each current under consideration are given by: 
Vector Current |4j] (T — 7 M ) 



2V, 



v 



m(3 — u 2 



III 



v 



ii 



1 - — A(u) + P v (u) In 



l + u 
1 - u 



+ Qv(u) 



P v (u) = — (33 + 22m 2 - 7m 4 ) 
Q v (u) = -{5u-3u 3 ). 



(B4) 
(B5) 

(B6) 

(B7) 
(B8) 
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Axial-vector Current @ (r = 1blu{q^qu/ q 2 ~ 9^)) 



N A 

K 



u 



l + u 



u 2 A(u) + P A (u) In j—^- + Q A 



u 



P A {u) = —(21 + 59m 2 - 19m 4 - 3m 6 ) 

Qa{u) = tt;(-21« + 30u 3 + 3m 5 ). 
Id 



A' Current gfl (r = d„j 5 ) 



N A , = 2 
1 



iL A' 



-M v 



2 
2 

3 
1 



1 + M 

m 2 A(u) + P A '(u) In + Q A ,(u) 

1 — u 



P A ,( U ) = — (13 + 28u 2 + 17u 4 - 2u e 
16 



<5a'(w) 
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-39m + 47m 3 + 6m 5 ) 



Scalar Current [g7J (r =1) 



N s = l 
s 2 



n 4 =2 



tiM'(u)+P s («) In 



1 +u 
1-u 



P 5 (u) = 4(3 + 34m 2 -13m 4 ) 
16 

Q s (u) = i(21w-3w 3 ). 



Qs(u) 



(B9) 
(BIO) 

(Bll) 
(B12) 
(B13) 

(B14) 
(B15) 

(B16) 

(B17) 

(B18) 

(B19) 
(B20) 

(B21) 

(B22) 

(B23) 



Pseudo-scalar Current [f47|l (T = 75) 

N P = 1 
3 



TT° 

Hp- 2 
2 



it 



n 4 P 



1-1-7/ 

A / (M) + P P (M)ln-^ + g P (M) 
1 — u 



Pp(u) 



1 

16 
1 



(19-48u + 2m 2 + 3u 4 



Qp{u) = ^(21u-3m 3 ) 

o 



(B24) 
(B25) 

(B26) 

(B27) 

(B28) 
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APPENDIX C: 



Resummation of Coulomb Singularities 

In this appendix we justify the resummation procedure used in Eq. ( |29| ) to avoid the 
adverse effects of threshold Coulomb singularities in the moment ratios for large n. 

As is well known [[HI] (see also the perturbative series for Imll(s) is generically an 



expansion in powers of A/ii, where u is the quark velocity introduced in Section 3. The series 
suffers from threshold singularities for small u, and this is accentuated in moment space, 
where the dominant contribution to the n th moment comes from the region u ~ 1/ ^Jn. In 
other words, for higher order moments the system appears more and more nonrelativistic, 
and the perturbative series is poorly behaved; the expansion parameter is now ^/n\. 

This behaviour is usually accounted for by treating the nonrelativistic Coulomb system 



exactly, leading to a resummation, e.g. in the form of the Sommerfeld-Sakharov factor [42 



ImU(X) — > " . with X = — -. (CI) 

v ; l-exp(-A) 3 u K J 

Since our conformal mapping procedure makes use of the same perturbative coefficients, a 
resummation of this form is still necessary. However, the modified virtuality distribution 
introduced in Section 2 actually allows an alternative resummation to that mentioned above. 

In order to illustrate the procedure, we consider a simple model and study first the way 
in which the standard perturbative approach breaks down for large n. Consider the n th 
moment, 



OC 



M "< r > = /„ d "j7+Yf^' < C2 > 

where a = s — 4m 2 is the virtuality, and Y = Q 2 + Am 2 . Now assume that the full spectral 
density p(a) is known and has the form 



**> = Ut^)- (C3) 

with p a positive integer, and normalized to p(a, A = 0) = 1. A perturbative expansion of 
p(a) has the form 

l-f- + o(-V, (C4) 

and thus exhibits Coulomb singularities for small u. 

For large n we may rewrite the moments in the formQ, 

M n (Y) ~ | o °° dap(a) exp (-n£) (l + oQ), (C5) 




8 We assume that p(a) is smooth 
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and observe that the dominant contribution to the integral corresponds to virtualities in the 
region a ~ Y/n (this is just the threshold behaviour u ~ noted earlier). Then, for 

the exact moments, M^ x , given by ( |C2|) with the exact spectral density (|C3|), the dominant 
threshold behaviour at small u (or large n) corresponds to 

p ex — > ( 1 + VZ-$=) , (C6) 



Y , 

which is positive and finite for all n. 

If instead we consider the 0(A) perturbative contribution only, we find 

Ppert — >l-pVn-^=, (C7) 
V Y 



which implies that the approximation must break down, and hence the moments will become 
unreliable, for large n when p pert passes through zero and becomes negative. 

We can now introduce a resummation which restores the correct behaviour of the mo- 
ments for large n. Performing the integration by parts which ensures consistency of the 
perturbative ambiguities with the OPE naturally suggests that we subsume higher order 
corrections to the spectral function into the denominator of the virtuality distribution. If 
we consider this as a formal transformation, we rewrite the moments in the form 



where the function f(cr') obeys the implicit equation 



rfW) 

a' = / dap(a). (C9) 



o 



This representation follows from considering the change of variables a = f(cr') and assuming 
that p(a) satisfies certain conditions which ensure the limits /(oo) = oo and /(O) = 0. 

Although this transformation is exact, in this paper we work to 0(A) and, solving (C9) 
to this order, we find 

f{a') =a' + 2p\V^. (CIO) 

Note that at this order, the present resummation essentially coincides with the Sommerfeld- 
Sakharov factor (|CT|). The corrections are of 0(X 2 /u 2 ) and are thus neglected. Inserting the 
leading order result in (|U8|) , and representing it in the standard form ( |C2| ) with a spectral 
density p re (a'), one finds that for small u ~ l/y/ri, 

Pre — > ( 1 + 2v^4=l • (Cll) 



Y t 

This clearly indicates an approximate recovery of the correct large n dependence associated 
with p ex (a), despite using only the 0(A) terms in p(a). 

To illustrate the success of this resummation we plot moment ratios R n = M n /M n _i using 
the exact (-Rjf ), perturbative (RP ert ), and resummed (i?^ e ) spectral functions, in Fig. 6. 
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0.2 1 1 1 1 L- 

10 



20 

n 



30 



FIG. 6. Moment ratios are plotted using the exact (R^), perturbative (it^ er *), and resummed 
{R r n) spectral functions. We use the parameters Y = 1, A = 0.1, and p = 1. Other choices lead to 
similar qualitative results. 

One clearly observes the close agreement between and as n becomes large, while 
the perturbative approximation without any resummation breaks down for large n. 

Furthermore, it is important to note that in the resummed case u ~ lii n °t the 

dominant contribution to the integral. This representation allows a saddle point to appear 
for a finite value of a which is independent of n. This is associated with the shift in the 
peak of the virtuality distribution observed in Fig. 1 for the massless case. Therefore, in 
practice the dominant contribution to the moments does not move closer to the threshold 
region, and as n becomes large the saddle point dominates the moments, leading to a finite 
limit. This behaviour is also observed in the QCD discussion presented in Section 3. 
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